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Abstract 

The van der Waals coefficients for the alkali-metal atoms from Na to Fr 
interacting in their ground states, are calculated using relativistic ab initio 
methods. The accuracy of the calculations is estimated by also evaluating 
atomic static electric dipole polarizabilities and coefficients for the interaction 
of the atoms with a perfectly conducting wall. The results are in excellent 
agreement with the latest data from ultra-cold collisions and from studies of 
magnetic field induced Feshbach resonances in Na and Rb. For Cs we provide 
critically needed data for ultra-cold collision studies. 
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The van der Waals interaction plays an important role in characterizing ultra-cold col- 
lisions between two ground state alkali-metal atoms. While the calculation of interaction 
coefficients has been a subject of great interest in atomic, molecular and chemical physics for 
a very long time, it is only very recently that novel cold collision experiments, photoassocia- 
tion spectroscopy, and analyses of magnetic field induced Feshbach resonances have yielded 
strict constraints on magnitudes of the coefficients. Moreover, due to the extreme sensitivity 
of elastic collisions to the long-range part of the potentials, knowledge of the van der Waals 
coefficients influences predictions of signs and magnitudes of scattering lengths. Although 
many theoretical methods have been developed over the years to calculate van der Waals 
coefficients, persistent discrepancies remain. 

In this paper, various relativistic ah initio methods are applied to determine the van der 
Waals coefficients for the alkali- metal dimers of Na to Fr [|1|] . As a check on our calculations, 
we also evaluate the atom-wall interaction constants, which have recently been calculated 
by other methods, and use them as a sensitive test of the quality of our wave functions. 
Furthermore, we calculate atomic polarizabilities and compare them to experimental data, 
where available. 

The dynamic polarizability at imaginary frequency aiiiS) for a valence state \v) can be 
represented as a sum over intermediate states \k) 

«m = §s , P ^"f^ M m ■ m\v) , (i) 



(E k - E v y + lu 



where the sum includes an integration over continuum states and R = J2j=i r j is the dipole 
operator for the iV-electron atomic system. We use atomic units throughout. The dispersion 
coefficient Cq of the van der Waals interaction between two identical atoms is 

C 6 = - / du Hiuj)} 2 . (2) 

71 Jo 

The coefficient C3 of the interaction between an atom and a perfectly conducting wall is 

1 f°° 

C3 = — / du a(iuj) , (3) 
An Jo 

or alternatively 

C 3 = ±(v\R-B.\v). (4) 

Using the latter relation, we have previously @ determined the values of C3 coefficients for 
alkali-metal atoms using many-body methods. 

The dipole operator R, being a one-particle operator, can have non-vanishing matrix 
elements for intermediate states represented by two types of Slater determinant. Firstly, 
the valence electron v can be promoted to some other valence state w. Secondly, one of the 
core orbitals a can be excited to a virtual state m, leaving the valence state v unchanged. 
In the language of second-quantization, the first type of states is represented by |0 C ) and 
the second type by a^aaO^ |0 C ), where |0 C ) describes the core. These states will be referred 
to as "valence" and "autoionizing" states, respectively. 
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In accordance with such a classification, we break the total polarizability a into three 
parts: the polarizability due to valence states a v , the core polarizability a c , and the valence- 
core coupling term a cv , with 

a = a v + a c + a cv . 

The last two terms arise from the summation over autoionizing states. In evaluating the 
core polarizability we permit excitations into all possible states outside the core. The term 
counter term accounting for the consequent violation of the Pauli principle. 
Various states contribute at drastically different levels to the dynamic polarizability. 
For example, 96% of the static polarizability of Cs is determined by the two intermediate 
valence states 6P1/2 and 6P3/2, other valence states contribute less than 1%. The core 
polarizability accounts for approximately 4% of the total value and the contribution of the 
core- valence coupling term is about —0.1%. The relative sizes of contributions to the static 
polarizabilities for the other alkali-metal atoms are similar. The dynamic polarizability 
a(iu>), given in Eq. ([[]), behaves as 



a{iu) ~ fvk/u 2 = N/u 2 , 



at large value of u, where we have used the nonrelativistic oscillator strength sum rule 
S (0) = fvk = N. Because the ratio a c J a v nonrelativistically is close to N — 1 we 
expect the core polarizability to give the major contribution at large 10. Therefore, the core 
polarizability becomes increasingly important for heavier atoms. 

Based on the above argument, we use several many-body techniques of varying accuracy 
to calculate the different contributions to the total polarizability. In particular, we employed 
the relativistic single-double (SD) all-order method to obtain the leading contribution from 
valence states ||. The core polarizability is obtained from the relativistic random-phase 
approximation (RRPA) ||. The core- valence coupling term and the non-leading contribu- 
tion from valence states is estimated in the Dirac-Hartree-Fock approximation by a direct 
summation over basis set functions ||. 

The relativistic single-double (SD) all-order method has been previously used to obtain 
high-precision atomic properties for the first few excited states in alkali-atom systems j3|. 
The results of theoretical SD matrix elements and comparison with experimental data are 
presented elsewhere ||. Generally, the electric-dipole matrix elements for principal transi- 
tions agree with precise experimental data to better than 0.5% for all alkali-metal atoms; 
the calculations being more accurate for lighter elements. In the present work, for Na, K, 
Rb, and Cs, we have used SD matrix elements for the first six lowest P\ii and P3/2 levels. 
For Fr, we have used SD matrix elements for a principal transition and matrix elements cal- 
culated with the third-order many-body perturbation theory (MBPT), described in JF], for 
four other lowest Pi/ 2 and P 3 / 2 states. Unless noted otherwise, we have used experimental 
values of energy levels from Ref. |8j and from the compilation of Dzuba et al. || for Fr. 

The relativistic random-phase approximation (RRPA) was used previously to obtain 
static core polarizabilities for all alkali-metal atoms except Fr in Ref. Q. In the present 
calculations we reformulated the original differential equation method used in H in terms 



of basis sets M, in a manner similar to [10|. We reproduce the results of Ref. [§J and, in 



addition, obtain a value of 20.41 a.u. for the static dipole polarizability of the Fr + ion. Zhou 
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and Norcross |Tl] find a c (0) = 15.644(5) for the polarizability of Cs + , by fitting Rydberg 
states energies to a model potential for Cs, while the present RRPA calculations yield the 
value a c (0)=15.81. Based on this comparison, we expect the RRPA method to give at least 
a few per cent accuracy in the calculation of a c (iw). 

To demonstrate the sensitivity of our results to errors in the core polarizability, we present 
the ratios of values calculated omitting a c to the total values of a(0), C3, and Cq in Table |. 
We see that while ct(0) is affected at the level of a few per cent, the core contribution to 
Cq becomes increasingly important for heavier systems. a c (iw) contributes 2% to Cq for 
Na and 23% for Fr. The atom-wall interaction constant C3, obtained with Eq. (|3|), is the 
most sensitive to the core contribution. Indeed, while a c contributes 16% of C3 for Na, it 
accounts for the half of the total value of C3 for Fr. 

The tabulation of our results for static dipole polarizabilities, atom-wall interaction con- 
stants C 3 , and Cq dispersion coefficients is presented in Tables PHPy]. In Method I we use 
high-precision experimental values [O for dipole matrix elements of the principal transi- 



tion. We used a weighted average of experimental data if there were several measurements 
for a particular transition. In Method II we use the theoretical SD matrix elements for the 
principal transition. We recommend using the values obtained with Method I for a(0) and 
Cq, since the accuracy of experimental data for the principal transitions is better than that 
of SD predictions. 

In Table || we compare our calculations with experimental data for static polarizabilities. 
We find perfect agreement with a high-precision value for Na obtained in recent atom- 



interferometry experiments []13|| . The experimental data for static polarizabilities of K, Rb, 
and Cs are known with the accuracy of about 2% [[TJ|,[nj. While we agree with those 
experimental values, we believe that our theoretical approach gives more accurate results, 
mainly due to the overwhelming contribution of the principal transition to the sum over 
intermediate states. The electric-dipole matrix elements for principal transitions are known 
typically at the level 0.1% accuracy for all alkalis. The theoretical error is estimated from 
the experimental accuracy of matrix elements [12], from an estimated 5% error for the core 
polarizabilities, and 10% error for the remaining contributions to a(0). 

A sensitive test of the quality of the present dynamic polarizability functions is obtained 
by calculating C3 coefficients in two different ways: i) by direct integration of a(iu) using 
Eq. (^) and ii) by calculating the diagonal expectation value of R 2 in Eq. (^). In the 
present work we extend calculations of the expectation value of R 2 || in the SD formalism 
to obtain C3 values for Rb, Cs, and Fr. In the Table [TIT], we compare the SD values for 
C3 with those obtained in p[ using MBPT. The difference of 7% for Cs and 10% for Fr 
between SD and MBPT values is not surprising, since the MBPT |7| underestimates the 
line-strength of principal transitions by a few per cent for Cs and Fr. To make a consistent 
comparison between the C3 values obtained by integrating a(iu) and by calculating the 
expectation value, we have used SD energies and matrix elements in Method II calculations 



in Table |TJ. These C 3 values agree to about 0.6% for Na, 1% for K and Rb, 2.5% for Cs, 
and 3.4% for Fr. At present, it appears no experimental data are available for comparison. 
We assume that most of the error is due to the RRPA method used to calculate the core 
polarizability. Therefore, the error estimates in Cq are based on the accuracy of experimental 



matrix elements for the principal transition [12| , and by scaling the error of core contribution 
from C 3 to Cq, using Table |. 
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The comparison of Cq coefficients with other calculations is presented in the Table [IV. 
For Na the results are in good agreement with a semi-empirical determination [16|. The inte- 
gration over a(iu) as in Eq. (0) has been most recently used by Marinescu, Sadeghpour, and 
Dalgarno |17| and by Patil and Tang In contrast to the present ab initio calculations, 
both works employed model potentials. In addition, Ref. [[r7|] used corrections to multipole 
operators to account for core polarization effects with parameters chosen to reproduce the 
experimental values of static polarizabilities, which for K, Rb, and Cs atoms are known from 
experimental measurements with an accuracy of approximately 2%. The major contribution 
in the integration of Eq. (|2|) arises from the region of to — and the integrand is quadratic 
in a(iu). Therefore, until more accurate experimental values for static polarizabilities are 
available, the predictions JlTj of Cq for K, Rb, and Cs have an inherent (experimental) ac- 
curacy of about 4%. Theoretical uncertainty of the method used in Ref. |T^] is determined, 
among other factors, by the omitted contribution from core polarizability as discussed in 
Patil and Tang gg 



Refs. 



used model-potential calculations with analytical repre- 
sentations of wave functions and with experimental energies. They used a direct summation 
method in Eq. ([!]). The contribution from the core polarizability was not included as can 
be seen from Eq. (3.4) of Ref. | 18fl . In fact, this formula in the limit of large us results in 
1/uj 2 instead of the correct limit ct(iu) — > N/uj 2 , which follows from the oscillator 



strength sum rule. Therefore, the model-potential calculations generally underestimate the 
Cq coefficients. Indeed, from the comparison in Table |IV|, one can see that the Cg values 
from Ref. 



17] and Ref. [18| are systematically lower than our values. 



Maeder and Kutzellnigg pOf used a method alternative to the integral Eq. (|j) to calcu- 
late dispersion coefficients by minimizing a Hylleraas functional providing a lower bound. 
However, their prediction depended on the quality of the solution of the Schrodinger equa- 
tion for the ground state. For alkali-metal atoms, model potentials were used to account for 
correlations. The predicted static polarizabilities are several per cent higher than experi- 
mental values, and are not within the experimental error limits. However, for Cg coefficients 
we generally find good agreement with the values of Maeder and Kutzellnigg . 

Recently Marinescu et al. [^l|] presented calculations of dispersion coefficients of different 
molecular symmetries for Fr, using a model potential method similar to Ref 
in Table IV our result for Fr is significantly larger than the result of Ref 
this may be because the method of Ref. [^1|] does not completely take into account the 
contribution of the core polarizability, which accounts for 23% of Cq for Fr. 

Elastic scattering experiments and photoassociation spectroscopy have sensitively con- 

m 



171 



. As shown 
We believe 



re- 



strained the possible values of Cg for Na and Rb. Van Abeelen and VerHaar 
viewed spectroscopic and cold-collision data for Na, including data from recent observa- 
tions of magnetic field induced Feshbach resonances 23[| . They considered values for Na of 
1539 < Cq < 1583 and concluded that Cg = 1539 gave the best consistency between data 
sets. Our result for Na using Method I is in particularly good agreement with this value. 
Photoassociation experiments |24]] for Rb limits the Cg coefficient to a range 4400-4900 a.u. 
and even more recently a study of a Feshbach resonance in elastic collisions of 85 Rb 
concluded Cq = 4700(50). Our value Cg = 4691(23) is in excellent agreement with this 
experiment. For Cs, knowledge of the value of Cq is critical for predictions of the sign of the 
elastic scattering length p6f, though it has been demonstrated the resulting cross sections 



are not particularly sensitive to the value of C 6 [E7|. For Fr, the paucity of other dimer data 
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constrains quantitative theoretical collisional studies for the near future. As photoassocia- 
tion experiments move beyond the alkali-metal atoms to other atoms with many electrons 
such as Sr [28J and Cr p9[] , it will be important to have reliable ab initio methods for cal- 
culation of atomic properties. The approaches presented here could, in principle, be applied 
to Sr and perhaps with some significant effort to Cr. 

AD would like to thank H. R. Sadeghpour, B. D. Esry, F. Masnou-Seeuws, and D. J. 
Heinzen for useful discussions. The work of AD, WRJ, and MSS was supported in part by 
NSF Grant No. PHY 95-13179 and that of JFB by NSF Grant No. PHY 97-24713. The 
Institute for Theoretical Atomic and Molecular Physics is supported by a grant from the 
NSF to the Smithsonian Institution and Harvard University. 
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TABLES 



TABLE I. Demonstration of the relative importance of the contribution of autoionizing states 
with increasing number of electrons N, where a", C3, and Cg represent values calculated disre- 
garding autoionizing states. 





Na 


K 


Rb 


Cs 


Fr 


a v {0)/a{0) 


0.99 


0.98 


0.97 


0.96 


0.94 


C$/C 3 


0.84 


0.73 


0.65 


0.59 


0.50 


Cq/C 6 


0.98 


0.93 


0.89 


0.85 


0.77 



TABLE II. Comparison of static dipole polarizabilities a(0) for alkali-metal atoms in atomic 
units. Method I designates the use of high-accuracy experimental data for electric-dipole matrix 
elements for principal transition. Method II designates the use of all-order SD values instead. 







Na 


K 


Rb 


Cs 


Fr 


Method I 
Method I 
Expt. 1 13 
Expt. [14 


a 
I 

15 


162.6(3) 

163.0 

162.7(8) 

I b 


290.2(8) 
289.1 

293.6(6.1) 


318.6(6) 
316.4 

319.9(6.1) 


399.9(1.9) 
401.5 

403.6(8.1) 


317.8(2.4) 
315.1 



a Values recommended from the present work. 

b Weighted average of experimental data from Refs. | p^Jl~5| . 



TABLE III. Comparison of atom-wall interaction constants C3 for alkali-metal atoms in atomic 
units. Method I designates the use of high-accuracy experimental data for electric-dipole matrix 
elements and energies for principal transition. Method II designates the use of all-order SD values 
instead. 





Na 


K 


Rb 


Cs 


Fr 


Method I, Eq. © 
Method II, Eq. (|) 
±(R 2 }, SD ab , Eq. @ 
i(JZ 2 ), MBPT §, Eq. (|) 


1.871 
1.875 
1.8858 
1.8895 


2.896 
2.887 
2.860 
2.838 


3.426 
3.410 
3.362 
3.281 


4.269 
4.247 
4.143 
3.863 


4.437 
4.427 
4.281 
3.870 



a Values recommended from the present work. 

b The values for Na and K are from Ref. ^], and those for Rb, Cs, and Fr are the present 
calculations. 
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TABLE IV. Tabulation and comparison of C% dispersion coefficients for alkali-metal atoms in 
atomic units. Method I designates the use of high-accuracy experimental data for electric-dipole 
matrix elements for principal transition. Method II designates the use of all-order SD values 
instead. 



Na 



K 



Rb 



Cs 



Fr 



Method I a 
Method II 



Ref. p| 



Ref. 120 



Ref. |7| 



Ref. |18| 



1556(4) 

1564 

1561 

1540 

1539 c 

1500 



3897(15) 
3867 

3945 
3813 
3796 



Ref. |1| 
Expt. p4| 



Expt. H 



4691(23) 
4628 

4768 
4426 
4531 

4400-4900 
4700(50) 



6851(74) 
6899 

6855 
6331 
6652 



5256(89) 
5174 



3934 d 



a Values recommended from the present work. 
b Semiempirical method. 

c For Na the value from Ref. |17| is 1472, obtained using the data from Ref. [14|. Using the same 



method, but with data from Ref. [13], the resulting value is 1539 [JTt 
d Value for 3 1 U molecular symmetry. Values for other symmetries are C$(Hf£ 
= 3947. 



3929, and C 6 (\) 
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